Effect of the lattice misfit on the equilibrium shape of strained islands in 

Volmer- Weber growth 



Jose Emilio PrietcQ 

Centro de Microandlisis de Materiales ( CMAM), 
Dpto. Fisica de la Materia Condensada and Instituto "Nicolas Cabrera", 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain 

Ivan MarkcnQ 

Institute of Physical Chemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria 

(Dated: October 15, 2010) 

We have studied the effect of the misfit on the equilibrium shape of three-dimensional pyramidal 
islands grown on a foreign substrate in the case of incomplete wetting (Volmer- Weber mode of 
growth). By means of atomistic simulations using anharmonic interaction potentials, we find that 
tensile islands have smaller aspect ratios compared with compressed islands owing to their better 
adhesion to the substrate. The average strains of consecutive layers decrease faster with thickness in 
compressed than in tensile islands. The strains decrease rapidly with thickness, with the consequence 
that above a certain height, the upper layers of the pyramid become practically unstrained and do 
not contribute to a further reduction in the top base. As a result, the truncated pyramids are not 
expected to transform into full pyramids. Our results are in good agreement with experimental 
observations in different systems. 
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Elastic strains play an important role in the growth 
of heteroepitaxial thin hims. In the case of complete 
wetting of the substrate by the film material, the elastic 
strain due to the lattice misfit results in a transition from 
planar to three-dimensional (3D) island growth beyond 
a critical thickness. i 2 - This is the well known Stranski- 
Krastanov (SK) mode of growth, which is widely used to 
produce quantum dots in semiconductor systems such as 
Ge/S: 3 and InAs/GaAs 4 . In the case of incomplete wet- 
ting (non-zero wetting angle), 3D islanding takes place 
from the very beginning of the deposition, irrespective 
of the value of the misfit. Nevertheless, in this case of 
Volmer- Weber (VW) growth, the misfit can play a crucial 
role in determining the equilibrium shape and, in partic- 
ular, the aspect ratio of the 3D islands. The two cases 
mentioned above differ considerably. In the SK mode, 
complete wetting and misfit strain operate in different di- 
rections: complete wetting favours planar growth, while 
strain favours clustering. This gives rise to an instabi- 
lity of planar growth against 3D islanding. In the VW 
mode both the incomplete wetting and the effect of misfit 
strain favor islanding. 

For the reasons given above, the equilibrium shape of 
islands on misfitting substrates has been the object of 
intensive studies^—. The dependence of the equilibrium 
shape in both SK and VW modes has been considered 
in detail by Muller and Kern* 1 ^— In the case of VW 
growth, they found that box-shaped crystals thicken with 
increasing value of the misfit. A qualitatively similar con- 
clusion was reached in Ref. 15[ under the condition that 
the crystal does not relax and is homogeneously strained 
to fit the substrate. In the case of pyramidal islands, 
the above authors found that the upper base shrinks and 
the island shape can change from truncated to complete 



pyramids with increasing misfit. It has been found in ad- 
dition that the equilibrium shape depends on the square 
of the lattice misfit thus being independent of the misfit 
sign. 14 Muller and Kern explained this result by neglect- 
ing surface stress effects in their calculations. 

In the present paper we study the effect of the lattice 
misfit, and in particular, of its sign on the equilibrium 
shape of pyramidal crystals. We performed atomistic 
simulations considering square-based pyramidal islands 
with fee atomic structure and (100) orientation, located 
on a substrate to which periodic boundary conditions are 
applied. The equilibrium shape of fee crystals consider- 
ing only nearest neighbor interactions is a truncated oc- 
tahedron or a cube-octahedron enclosed by six (001) and 
eight (111) faces with equal edge lenghtsJ^ The wetting 
function (or adhesion parameter) is defined a a 17 ' 18 



(1) 



where ip' and ip are the substrate-deposit (adhesion) and 
deposit-deposit (cohesion) atomic interaction energies, 
respectively. If the wetting function is smaller than 0.25 
at zero misfit, the lateral cubic faces do not appear and 
the equilibrium shape remains a simple truncated pyram- 
id (Fig.[T]). The lattice misfit decreases the degree of wet- 
ting and contributes to an effective increase in the wet- 
ting function^^ For this reason, we restrict ourselves 
to values of the wetting in the interval <j> = -j- 0.2. 

We have performed simple minimization calculations. 
The atoms in the islands interact through a pair poten- 
tial whose anharmonicity can be varied by adjusting two 
constants fi and v (/z > v) ^ > 22 
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where b is the equilibrium atom separation. For ^ = 2v 
the potential adopts the familiar Morse form. It is worth 
to note that these kinds of potentials, due to their spher- 
ical symmetry, are better suited to describe the bonding 
in metals. By contrary, semiconductors are characterized 
by directional, covalent bonds and bonding angles need 
to be taken into account. However, the potentials ap- 
propriate for semiconductors are also anharmonio 2 ^, so 
we expect our results using anharmonic potentials to be 
qualitatively valid for those materials as well. 

Energy minimization is performed by allowing the 
atoms to relax until the forces fall below some negligible 
cutoff value. In spite of its simplicity, the above potential 
includes the necessary features to describe real materials 
(theoretical strength and anharmonicity) by varying the 
constants /i, v and Vq. We consider interaction only in 
the first coordination sphere and simulate both rigid and 
relaxable substrates, in which a given number of layers S 
are allowed to relax. The initial positions of the atoms 
correspond to the centers of the potential troughs pro- 
vided by their neighbors underneath; those neighbors are 
not updated when atomic displacements approach val- 
ues close to one half of the interatomic distance, so our 
model is not appropriate for the description of configu- 
rations containing misfit dislocations. The interatomic 
spacing in the substrate material is a so that the lattice 
misfit is given by / = (b — a) /a. The case of positive mis- 
fits corresponds to the lattice parameter of the deposit 
being larger than that of the substrate, so the material 
in the islands is compressed. The atoms interact with 
their lateral neighbors mainly through the stronger re- 
pulsive branch of the potential. In the opposite case, the 
islands are tensile and the atoms interact with their lat- 
eral neighbors with the weaker attractive branch of the 
potential. As a result, the surface stress of the side walls 
should differ in compressed and tensile crystallites, or in 
other words, they should relax in different degrees. The 
side walls of compressed islands are expected to relax by 
a greater amount compared with tensile islands. 

In addition, the misfit sign affects the adhesion of the 
crystals to the substrate in a different way. In compressed 
islands the atoms of the first atomic layer are strongly 
displaced from the bottom of the corresponding potential 
troughs of the substrate owing to the stronger repulsion 
of the lateral neighbors. As a result, the adhesion is 
reduced. On the contrary, the first-layer atoms in tensile 
islands are only slightly displaced from the bottoms of 
the respective potential troughs of the substrate and the 
adhesion is stronger compared with compressed islands 
at the same absolute value of the lattice misfi t 26 ' 27 . 

In order to determine the equilibrium shape of the is- 
lands, we apply the following procedure. We consider is- 
lands of different heights and bases with shapes close to 
squares, either of n x n or of n x (n + 1) atoms, where n is 
an integer. Allowing the structures to relax, we calculate 
the total energy per atom of the islands, for heights h 
varying by 1 monolayer (ML), as a function of the total 
number of atoms N in the pyramid. Due to the dis- 




FIG. 1: (Color online) Epitaxial fcc(lOO) island of -7% misfit 
on a rigid substrate. It has the shape of a truncated square 
pyramid with 30 x 30 atoms in the base plane and a height of 
4 monolayers. The color scale denotes the level of hydrostatic 
strain as given by the atomic local strain tensor— and has 
been represented using the AtomEye software—. Strain is 
highest at the center of the island and at the bottom layer 
and it is relaxed at island edges and at the topmost layer. 



Crete character of the island sizes produced in this way, 
we need to interpolate between the two closest values at 
the corresponding height to the desired total number of 
atoms, e.g. N = 1000. In this way, we find the height of 
the island with lowest energy. For example, Fig. [5] shows 
that four-layers thick islands have lowest energy when the 
number of atoms is smaller than 300. At N > 300 the 
lowest energy thickness is 5 MLs. We then calculate the 
lowest energy heights of islands at a given constant num- 
ber of atoms and for different values of the lattice misfit. 
Once we know the lowest energy heights we determine 
the edge lengths of the lower, n, and upper, n', bases 
from the relation n' = n — h+l and from the formula for 
the total number of atoms N(n, h) = Y^k=i ^ 2 ~ Y^k=i ^ 
as a function of n and h. The height h is measured in 
number of monolayers and n and n' give the number of 
atoms in the upper and lower edges, respectively. Thus, 
the full pyramid is characterized by n' = 1 and h/n = 1. 

The simplest expression for the equilibrium shape of 
the pyramid is h/n' = 40, where (f> is the wetting func- 
tion defined above and the coefficient 4 reflects the coor- 
dination of an atom on the (001) surface. This formula 
has been derived by comparing the works per atom to 
evaporate whole atomic planes of the lateral pyramidal 
faces and the upper base, a method which had been in- 
troduced by Stranski and Kaischew 2 ^ (an excellent recent 
review is given in Ref. [29t ). However, we will illustrate 
our results in terms of plots of the ratio h/n vs. the 
lattice misfit because the ratio h/n' tends to the island 
height h, which depends on the number N of atoms in 
the pyramid when it becomes complete and the upper 
base disappears, while the ratio h/n tends to unity. 

Figure [3] demonstrates our main result. It shows the 
dependence of the aspect ratio h/n of islands of mini- 
mum energy on the value of the misfit, both negative and 
positive ones. The islands consist of small (a) N = 250 
and large (b) N = 1000 number of atoms. The substrate 
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FIG. 2: (Color online) Dependence of the total energy per 
atom as a function of the number of atoms TV in the pyra- 
mid for islands with different heights varying through one 
monolayer. As seen islands with TV < 300 have lowest energy 
when they are four monolayers high whereas the lowest energy 
height is 5 monolayers at TV > 300. The wetting parameter is 
<j> = 0.2, the lattice misfit / is zero and a rigid substrate was 
assumed in the simulations. 



is assumed rigid. Several important conclusions can be 
reached. First, the lattice misfit has a great effect on 
the equilibrium shape. The aspect ratio increases steeply 
beyond a moderate misfit of about 3%. The smaller the 
wetting function (i.e. the stronger the adhesion), the 
stronger the effect of the misfit. Or, the stronger the 
tendency to 3D islanding due to the weaker wetting, the 
smaller is the effect of the lattice misfit. The effect of 
the misfit is greatest at <fr — >• 0. As mentioned above, the 
reason is that both effects are complementary. 

Second, it is seen that a positive misfit has a percep- 
tibly greater effect on the equilibrium shape than a neg- 
ative one of the same absolute value. This is due to the 
stronger repulsive forces of the interatomic potential and 
to the weaker adhesion. In tensile overgrowth the in- 
crease in the aspect ratio is delayed to greater absolute 
values of the misfit. And more important, increasing the 
lattice misfit does not lead to disappearance of the upper 
base, h/n — > 1, as predicted by Miiller and Kerr>i4 In 
both cases of positive and negative misfit the aspect ratio 
h/n goes to a saturation value of about 0.75 irrespective 
of the island size and misfit sign. 

The curves h/n vs misfit show an initial increase and 
an inflection point. Beyond this point the curves tend to 
a saturation value lower than unity as mentioned above. 
Note that at large absolute values of the misfit the in- 
terface between the crystallites and the substrate is ex- 
pected to be resolved in a cross grid of misfit dislocations. 
The arrows in Fig. [3] show the approximate values of the 
misfit at which dislocations are introduced in compressed 
islands. The tensile islands, on the other hand remain 
coherent. Although our model, as explained above, is 
not well suited for the description of dislocated config- 
urations, it seems clear that the tendency to saturation 
begins long before the introduction of dislocations. Fur- 
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FIG. 3: (Color online) Misfit dependence of the aspect ratio 
h/n for different values of the wetting function 4> = 0.0, 0.03 
and 0.10 and a total number of atoms in the pyramid of (a) 
TV = 250, (b) TV = 1000. A rigid substrate was assumed in 
the calculations. The arrows mark the appearance of misfit 
dislocations in the growing islands. 



thermore, since both 3D clustering as well as the intro- 
duction of misfit dislocations are (possible competing) 
mechanisms that contribute to the relaxation of epitax- 
ial strain, we can conclude that the introduction of misfit 
dislocations is not the basic reason for the islands to re- 
main truncated pyramids at the sizes studied so far. As 
will be shown below the main reason is the fast decrease 
in the strain with the pyramid height. 

We can understand the different behavior of the as- 
pect ratio at positive and negative misfits as follows. As 
shown in 1+1 dimensional model s 19 ' 30 , the thermody- 
namic driving force for 3D islanding in coherent Stranski- 
Krastanov growth is the reduced adhesion of the initial 
2D islands to the wetting layer. This is due to the dis- 
placements of the atoms located closer to the islands 
edges from the respective potential troughs of the sub- 
strate (the wetting layer). Thus the average adhesion 
parameter is non-zero and 3D islanding is thermodynam- 
ically favored. A more detailed study of the adhesion 
parameter as a function of the misfit sign has shown 
that compressed islands display larger values of the adhe- 
sion parameter (weaker wetting) compared with tensile 
islands. Hence, greater absolute values of the (negative) 
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misfit in tensile islands are required in order to reach the 
same values of the adhesion parameter as compared to 
compressed islands. 
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FIG. 4: Misfit dependence of the adhesion parameter, <j>'(f), 
for compressed and tensile islands for the case of (f> — 0.0 (SK 
mode of growth). Results are shown for different numbers S 
of substrate layers allowed to relax: circles for S = (i.e., 
rigid substrate) and squares for S = 3 ML. The islands are 
5 ML high and contain a total amount of N = 990 atoms. 

The same result is obtained in our 2+1 dimensional 
model as demonstrated in Fig. @] for the Stranski- 
Krastanov growth. Here the wetting parameter <f> due 
to different bonding, as defined by Eq.([T]), equals zero 
and we can define a mean wetting parameter <p' as 
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where t/ a( jh is the mean interaction energy between an 
atom at the bottom layer of the island and the substrate, 
which equals -4Vb at zero misfit. The tensile islands al- 
ways display better adhesion to the substrate. Although 
the difference in adhesion is greater in rigid substrates 
compared with relaxed ones, the effect of the potential 
anharmonicity remains significant. The adhesion para- 
meter for tensile islands remains close to zero (complete 
wetting). 

In the case of Volmer- Weber growth the effective ad- 
hesion parameter in misfitting overgrowth, c g, is thus 
composed of two parts. The first part is due to the differ- 
ence in bonding and is given by Eq. ([lj , while the second 
part is due to the lattice misfit as discussed above. We 
can thus write 



<M/) = 4> + <t>'{f) 
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where the second term in the right-hand side depends 
solely on the lattice misfit, /, so that <fi'(0) = 0. The 
inequality <j>'(+\f\) > </>'(— and in turn e ff(+|/|) > 
0eff (— explains the results shown in Fig. [3J 

We now study the strain distribution in tensile and 
compressed islands on a substrate which is allowed to 
relax. Fig. [SJa) shows the distribution of the strains in a 
compressed crystallite (/ = +4%). The crystal is 5 MLs 



thick and the same number of monolayers of the substrate 
are allowed to relax. As seen, the 4th and the 5th MLs 
of the substrate are very weakly strained. The strains in 
the substrate below the island do not exceed 1% whereas 
those in the first layer of the crystallite are close to -3%. 
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FIG. 5: (Color online) Distribution of the lateral bond strains 
in the center of the crystal and in the layers of the substrate 
underneath for misfits of (a) f = 4% and (b) f = -4%. Strains 
are referred to the lattice parameters of crystal a and sub- 
strate b, respectively. The wetting parameter due to the dif- 
ferent chemical nature of substrate and adsorbate is <f> = 0.03, 
islands of 5 ML height and 990 atoms were considered and 
5 ML of the substrate were allowed to relax. 

Figure E^b) shows the corresponding strain distribu- 
tion in tensile island at the same absolute value of the 
misfit (/ = —4%). The distribution is qualitatively the 
same. The strains in the island are again much larger 
than those in the substrate. This is easy to understand 
having in mind that the substrate crystal has no free 
edges and the substrate atoms are highly constrained 
by their lateral neighbors. As a result, the total strain 
energy stored in both compressed and tensile islands is 
about two orders of magnitude larger than the strain en- 
ergy stored in the substrate. The important conclusion 
we can extract is that the substrate behaves as relatively 
"stiff" against the epitaxial strain induced by the over- 
growth and can be safely considered rigid in comparison 
with the deformations of the islands. This justifies the 
results in Fig. [3] where the substrate was assumed to be 
rigid. 

A more careful inspection of Fig. [5] shows that the lat- 
eral strains of the separate layers in compressed islands 
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are always smaller than the corresponding strains in ten- 
sile islands. In addition, the strains of the separate lay- 
ers decrease with layer height, this decrease being larger 
in compressed than in tensile islands (Fig. [5]) . In other 
words, compressed islands tend faster to achieve the bulk 
lattice parameter of the overgrowth material compared 
with tensile islands. 

However, the height distribution of the strain energies 
of the separate layers shows an unexpected behavior (Fig. 
[7]). The strain energies of the layers beyond the second 
one in tensile islands become larger than those in com- 
pressed layers. This can be understood as follows. The 
more strained the layer the more significant is the effect 
of the potential anharmonicity. The first layers are highly 
strained and the steeper repulsive branch of the poten- 
tial leads to a very high strain energy of the compressed 
island. The tendency of atoms to relax the strain is thus 
higher in compressed overlayers and this leads to smaller 
values of the strain in spite of the higher amount of strain 
energy as compared to tensile islands at the same abso- 
lute value of the lattice misfit. Upper layers become less 
and less strained and the anharmonicity cannot overcom- 
pensate the smaller absolute values of the strains. As a 
result, the more strained upper layers in tensile islands 
possess also a larger amount of strain energy. 



they are expected to play a secondary role in preventing 
the transformation of the small islands into full pyra- 
mids. An additional argument is that misfit dislocations 
are introduced in tensile islands always at misfit values 
where saturation takes place. In order to further clarify 
the role of the anharmonicity of the potential, we also 
studied the island's aspect ratio as a function of the lat- 
tice misfit using harmonic potentials of the same spring 
constants as the anharmonic ones previously discussed. 
The curves (not shown) obtained for positive and neg- 
ative misfits nearly overlap and are located as expected 
between the curves corresponding to the anharmonic po- 
tentials, slightly closer to the curve of positive misfit. 
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FIG. 6: Average strains of consecutive layers in compressed 
and tensile pyramids with +4.0% and -4.0% lattice misfit, 
respectively. Islands are 5 ML high and contain 990 atoms, 
the wetting parameter is (j> = 0.03 and 5 ML of the substrate 
were allowed to relax. 

Figure [5] also shows that the first layer is practically 
pseudomorphic with the substrate (the strains are nearly 
equal to the misfit taken with opposite sign) with the 
exception of a few edge bonds. On the contrary, most 
of the bonds in the uppermost (fifth) layer are much 
less strained than the nominal misfit showing that the 
strains rapidly goes to zero. It follows that the upper 
layers of a complete pyramid are practically fully relaxed 
and so they do not contribute to the accumulation of 
strain in the crystallite. For this reason, the crystallites 
do not transform to full pyramids for reasonable values of 
the misfit as expected from the continuous and harmonic 
theory of Miiller and Kern. The presence of misfit dislo- 
cations favors further relaxation of the upper layes, but 



FIG. 7: Strain energy of consecutive layers in compressed 
and tensile pyramids with +4.0% and -4.0% lattice misfit, 
respectively. Islands are 5 ML high and contain 990 atoms, 
the wetting parameter is <f) — 0.03 and 3 ML of the substrate 
were allowed to relax. 

It is worth to note that the equilibrium aspect ratio of 
the 3D islands in the case of complete wetting (SK mode) 
follows the same trend as in the case of VW growth mode 
(Fig. 13]). The only difference is that the initial islands 
are single-layer islands (4> = 0) which then transform 
into multilayer islands^. The only physical reason for 
incomplete wetting is the lattice misfit which leads to 
displacements of the overlayer atoms from the bottom of 
the substrate potential trough o 19 ' 20 and thus to a weak- 
ening of the interfacial bonding. 

Experimental observations of equilibrium shape of is- 
lands on misfitting substrates, in the appropriate ranges 
of parameters to compare with the results of our model, 
are scarce in the literature. We will not discuss here the 
results concerning the equilibrium shape of metal clusters 
on MgO (for a review see Ref. [l6l ). One of the reasons 
is that wetting functions are greater than 0.5 and so the 
equilibrium shapes are not simple pyramids. Another 
reason is that all metals studied (Ni, Pd, Pt) show a neg- 
ative misfit with the substrate^. Nevertheless, smallest 
particles of Pd on MgO (f = -6%) studied showed an 
equilibrium shape of a truncated pyramid, h/n » 0.7—. 

Goldfarb et al*2i found that the equilibrium shape 
of small crystals of TiSi2 on Si(lll) is a flat hexag- 
onal pyramid with the (013) atomic plane parallel to 
Si(lll). The misfit in the two directions < 112 > and 
< 110 > amounts to 5% and 10%, respectively. Silly 
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and Castell studied a completely different system: Fe on 
SrTiO 3 (001) 35 . After annealing, the equilibrium shape 
of bcc-Fe islands is a square pyramid with an aspect ra- 
tio l/h = 1.2 irrespective of the island size, where I is the 
edge length of the upper base. The side walls make an 
angle of 45° with the substrate so that the aspect ratio 
h/L amounts to 0.31, where L is the edge length of the 
lower base. Nowicki et al. found that the equilibrium 
shape of Pb islands on Ru(001) corresponds to a trun- 
cated pyramid with the (111) face as the upper bas o 36 ' 37 . 
So in all cases discussed above, the equilibrium shape is a 
truncated pyramid in accordance with the results of the 
present paper. 

In conclusion, using anharmonic interaction potentials, 
we have found that compressed islands show larger val- 
ues of the aspect ratio than tensile ones at the same ab- 
solute value of the misfit owing to its larger effective ad- 
hesion parameter. Atoms belonging to the first layer of 
compressed islands show larger displacements from the 
bottom of the potential troughs of the substrate due to 
the stronger interatomic repulsive forces, which leads to 
weaker adhesion of the islands to the substrate. On the 



contrary, the weaker attractive branch of the potential 
acts between the atoms in tensile islands and the effec- 
tive wetting parameter is smaller. Due to the same rea- 
sons, compressed islands are less strained compared with 
tensile ones and tend to attend the bulk lattice param- 
eter more rapidly with increasing island thickness. As a 
result, the upper layers of the pyramids are practically 
unstrained and they cannot transform to full pyramids 
with increasing values of the lattice misfit. The strain 
energy stored in the substrate is about two orders of mag- 
nitude smaller than in the islands due to the constraint 
of the atoms by their lateral neighbors. As a result, for 
many calculations, the substrate can usually by assumed 
as rigid. 
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